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The ion toroidal rotation in a tokamak consists of an E x B flow due to the radial electric field 
and a diamagnetic flow due to the radial pressure gradient. The turbulent pinch of toroidal angular 
momentum due to the Coriolis force studied in previous work is only applicable to the E x B flow. 
In this Letter, the momentum pinch for the rotation generated by the radial pressure gradient is 
calculated and is compared with the Coriolis pinch. This distinction is important for subsonic flows 
or the flow in the pedestal where the two types of flows are similar in size and opposite in direction. 
In the edge, the different pinches due to the opposite rotations can result in intrinsic momentum 
transport that gives significant rotation peaking. 
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Introduction. Ion toroidal rotation can stabilize MHD 
resistive wall modes P, Q and increases the confinement 
time in tokamak plasmas Consequently, external 

momentum sources are frequently used to drive rotation 
in current tokamaks. However, for larger volume, high 
density plasmas such as those expected in ITER, the 
available external sources are insufficient to produce sig- 
nificant toroidal rotation. Thus, we can only rely on rota- 
tion that is "intrinsically" generated by the radial trans- 
port of the toroidal angular momentum due to microtur- 
bulence and by coupling to the rotation at the plasma- 
wall boundary. Radial turbulent transport redistributes 
the toroidal momentum within a tokamak. 

The intrinsic ion toroidal rotation has been measured 
in many tokamaks Q , and the size of the flow is much be- 
low the ion thermal velocity (Mach number ~ 0.1 — 0.2). 
The subsonic level of the flow is in accordance with the- 
oretical estimations 0, HJ . Due to the symmetry proper- 
ties of the turbulence in an up-down symmetric tokamak, 
the rotation intrinsically generated by momentum redis- 
tribution is subsonic when the rotation at the edge is sub- 
sonic, as it usually is. The reason is that for sonic flow a 
symmetry of the turbulence that gives vanishing toroidal 
angular momentum transport, and hence momentum re- 
distribution, can only be broken by pre-existing nonzero 
toroidal flow or radial flow shear of the flow, giving tur- 
bulent momentum pinch and diffusion, respectively 

The momentum pinch is inward momentum flux that 
depends on the flow size 13, 11|- For small flows, this 
momentum flux is simply proportional to the flow. It 
is not just convective transport due to the radial parti- 
cle pinch carrying momentum, but also momentum flux 
caused by the presence of non-zero rotation even in the 
absence of the particle pinch. 

In previous work about the momentum pinch 
12j, there is no distinction between the different rota- 
tion types. Both the radial pressure gradient and the 
radial electric field contribute to the subsonic toroidal 



rotation. The angular frequency of the toroidal rotation 
i\ has two pieces: one driven by the radial electric field, 
fl^E = — {c/RBg) (dd>o/dr), which corresponds to Ex B 
plasma flow; and the other driven by the pressure gra- 
dient, fi^p = — (c/ ZeriiRBg) (dpi/dr), corresponding to 
the diamagnetic flow. Here, c is the speed of light, R 
is the major radius, Be is the poloidal magnetic field, 
4>a{r) is the long wavelength electrostatic potential, r is 
the radial coordinate, Ze is the ion charge, rij(r) is the 
ion density, and Pi(r) is the ion pressure. The turbulent 
momentum pinch due to the radial electric field driven 
rotation is the only momentum pinch that has been in- 
vestigated by other authors [l2| . In this Letter, the tur- 
bulent momentum pinch for the toroidal rotation driven 
by the radial pressure gradient, i\. p , is calculated, and 
it is compared with the pinch for the rotation driven by 
the radial electric field, O^s. 

The two types of toroidal rotation have different ori- 
gins. The radial electric field is related to the toroidal 
rotation via the Lorcntz force, and as a result the ra- 
dial electric field changes on the momentum confinement 
time scale [l3| • Conversely, the pressure gradient is deter- 
mined by the turbulent transport of energy and particles, 
changing on the energy confinement time scale. 

We found that the effect of the two different types of 
rotation on the turbulent momentum transport are also 
different. The main difference is that the energy of the 
particles changes due to the radial electric field but not 
due to the radial pressure gradient. This is because the 
orbits of the particles in a tokamak deviate radially from 
the flux surface in which the magnetic field lies due to 
their magnetic drifts. The work on the particle done by 
the radial electric field is poloidally antisymmetric for an 
up-down symmetric tokamak, and it can give a nonzero 
momentum flux because it breaks the up-down symmetry 
of the turbulence. However, the pressure gradient cannot 
change the particle orbit or energy. 

The dependence of the turbulent momentum pinch on 
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the rotation type reveals a new symmetry breaking mech- 
anism in the absence of flow in which the two different 
types of rotation are equal in size and opposite in sign. 
This intrinsic rotation generation mechanism can be im- 
portant in the pedestal, a region of strong gradients ob- 
served in the edge of tokamaks in certain regimes 14j . 



where both a large pressure gradient and a strong ra- 
dial electric field contribute to the toroidal rotation in 
opposite directions. 

Radial momentum flux. If there is no external momen- 
tum source, the radial transport of ion toroidal angular 
momentum, II, determines the evolution of the toroidal 
flow, 
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where rrii is the ion mass, J is the Jacobian and (...) s is 
the flux surface average. In steady state, the radial mo- 
mentum flux has to be zero at every flux surface. The bal- 
ance between non-zero pieces in II gives the radial profile 
of toroidal flow; e.g., 11(17^, dVl^/dr, ...) = determines 
i\(r). The radial momentum flux is dominantly caused 
by the E x B radial drift due to the short wavelength 
fluctuating potential carrying the fluctuating toroidal an- 
gular momentum, 
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where Q and f are the toroidal and radial unit vectors, re- 
spectively, B and b are the magnitude and the unit vector 
of the magnetic field, respectively, y = b x f , and (...)t 
is a time average over several turbulence characteristic 
times. Here, v^ xB is the E x B drift due to the fluc- 
tuating electrostatic potential and V* h ^ is the fluctuating 
toroidal velocity. The fluctuating distribution function 
ftb and the fluctuating piece of the potential <j> tb due to 
microturbulence can be obtained from gyrokinetics. 

Gyrokinetics. Microturbulence in a tokamak is well de- 
scribed by the gyrokinetic approximation. Gyrokinetics 
assumes that the gyromotion is much faster than the tur- 
bulent fluctuations (u> <C f2j), and the turbulence char- 
acteristic length scales in the directions parallel and per- 
pendicular to the static magnetic field are of the order of 
the ion Larmor radius and the size of the device, respec- 
tively (p%/l± ~ 1 and lj_/l\\ ~ Pi/ a = P* *C 1), where 
lj is the frequency of turbulence, VLi is the ion Larmor 
frequency, pi is the ion Larmor radius, a is the minor 
radius, and In and l± are the parallel and perpendicu- 
lar turbulence length scale, respectively [15l [l6|. The 
dependence of the distribution function on the angle of 
Larmor gyration is eliminated by averaging over the fast 
Larmor gyration. As a result, the velocity space is de- 
scribed by two variables (the kinetic energy E = v 2 /2 



divided into two 
v-M\C|74 



and magnetic moment p — Vj_/B, where v± — ^Jv 2 — v 2 

and are the velocity perpendicular and parallel to the 
static magnetic field, respectively). After averaging over 
the gyromotion, the particle motion can be described by 
following the center of the Larmor gyration. 

We assume that the ion distribution function 
fi = fo L) + ftb ] can be 
pieces: — nij [tv'I 2 v\^ exp ( 

/m (l + ( mv \\/Ti) /B)) i s the lowest order non- 
fluctuating piece that includes the toroidal flow, and 
is the fluctuating short wavelength piece whose size is 
0(/5*/m)' Here, Jm is the Maxwellian distribution func- 
tion, / = B^R, Bq is the toroidal magnetic field and 
v t i = ^j2Ti/nii is the ion thermal velocity for the tem- 
perature Ti. We assume Bg/B -C 1 and that the toroidal 
flow is sufficiently small to make (B / Bg)p+vu ~ S\ iP i? ~ 
Q(,eR *C Da. Similarly, the electrostatic potential can 
be divided into two pieces (f> = <f>o + <t> th ■ The gyrokinetic 
equation for the turbulent piece of the ion distribution 
function in the lab frame in the presence of the ion rota- 
tion due to both radial electric field and pressure gradient 
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Here, v M = (p/ttijbxVB+ivpQi) 6x(6-V6) is the VB 

and the curvature drift, (...) is the average over the gyro- 
motion, and C is the ion collision operator. Equation ^ 
only includes the terms related to the diamagnetic flow 
among the higher order corrections of order (B / Bg)p 2 . 
Other terms of order (B/Bg) pi will be treated in a fu- 
ture publication 

EE 

1 1 71 1 . Notice that the acceleration 
and the drift of a single particle are not affected by the 
pressure gradient but are affected by the radial electric 
field — V<^0) while the background distribution function 
is modified by both the pressure gradient and the radial 
electric field. The gyrokinetic equations can be derived 
for every species and they are coupled by imposing the 
quasineutrality condition. 

The gyrokinetic equation with nonzero rotation can be 
described either in a lab frame (L) or a rotating frame 
(R). These frames are related by a rotating frame trans- 
formation, v'n = v\\ — IHq/B and e = E — IVl^v\\/B 13|, 
where v',, and e are the parallel velocity and the kinetic en- 
ergy variable in the rotating frame. In the rotating frame, 
the ion distribution function is fi = /q (v',> ) + f t ^ M, ) 
where the lowest order non-fluctuating piece is related to 



that in the lab frame by fQ R \v'^) = /m(vJ) = /q W («||), 
and the first order fluctuating piece is given by f^ b (uf, ) = 



3 



f { tb L \v\\) * /4 L) («||) + WB) (dflP/dv 1 ^). Applying 
this transformation to Eq. ((4]), the gyrokinetic equation 
in the rotating frame is obtained, 
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where = (p/ik) b x VB + (v'^/Q.^ b x (b ■ V6) is the 
V-B and the curvature drift in the rotating frame, and 
v' c = r2u[|0 f /QiJ 6 x [(V-R x C) x S] is the drift due to 
the Coriolis force. 

In the frame of rotation driven only by the radial elec- 
tric field (fi^B 7^ and f\ jP = 0), the gyrokinetic equa- 
tion is modified only by the additional Coriolis terms 
[l2j . without any energy derivative of the turbulent dis- 
tribution function. Conversely, the rotation driven by 
the pressure gradient (1\,.e = and f2f )P ^ 0) does not 
include the energy derivative terms in the lab frame but 
it does have these terms in the rotating frame. This 
additional acceleration term that depends on the type 
of the rotation, (ZeRBe/niic) v' M ■ Vr (dftb/de) re- 
sults in different momentum pinches for the two rotations 
in either lab frame Q or rotating frame ([5]). The accel- 
eration term breaks the symmetry of the turbulence by 
introducing an up-down asymmetry effect in the poloidal 
coordinate (0). For example, for a tokamak with circu- 
lar flux surfaces, the term is proportional to sin (9 which 
is odd in 0. The other rotation type dependent term, 
(RBg/B) Vr x b ■ V/t&fi^p, only gives a Doppler shift to 
the fluctuation that cannot change the radial transport. 

Numerical results. To model the momentum transport 
II(f2£, dft^/dr), we assume that is sufficiently small 
and linearize around = and dCl^/dr — 0, giving 
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where IT„ t is the intrinsic toroidal angular momentum 
flux in the absence of rotation and rotation shear (fi^ = 
and dSlt^/dr = 0), P^ is the pinch coefficient, and xc 1S 
the toroidal momentum diffusivity. We use the gyroki- 
netic code GS2 [l8| to evaluate the diffusion and pinch 
coefficients for the different types of rotation. In this Let- 
ter, we assume that the momentum diffusivity coefficients 
are not different for the different type of rotations and 
X(; is inferred from a fixed Prandtl number, the ratio of 
the momentum diffusivity to the ion heat flux diffusivity, 
(Pr = xc/x* = - (Xc/Qi) (nidTi/dr) =s 0.5), obtained 
in a simulation with velocity shear. 

The numerical results for the momentum pinches are 
given in Fig. 1. The default plasma parameters for the 
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FIG. 1: Time averaged ratio of ion toroidal angular momen- 
tum flux (II) to ion heat flux (Qi) as a function of rotation 
(£l( x Ro/vu) for zero rotation shear and the different types of 
rotation: radial electric field driven rotation (blue-star), pres- 
sure gradient driven rotation (green-circle) and opposite ro- 
tations of the two types (red-diamond). The error bars show 
the standard deviation of the fluxes from the time average 
values due to the typical turbulent fluctuations. 



simulations are Rq/Lt — 9.0, Rq/L u = 9.0, q = 2.5, 
r/a — 0.8, Ro/a = 3.0, and s — 0.8 where Lt 

/', ; '77 , fir] and L n = —rii/{dni/dr) are character- 
istic lengths of temperature and density, respectively, 
Rq ~ {R)s, q is the safety factor, and s is the magnetic 
shear. We use the rotation peaking factor Pq/xc to quan- 
tify the strength of the pinch (P(/xt = — (dQ^/dr) 
is the rotation peaking for II = and IT nt = in Eq. 
([6])). Fig. 1 shows that the rotation peaking factors are 

P c ,e/xc - 2.9/i? for H ( , E , and P (lP /xc ^ 3.5/ Ro for 
The difference is about 22% of the pinch. 

Interestingly, the difference of the pinches results in 
negative momentum flux even for zero rotation (f2^ = 
f2f,E + &C,p — 0, and p > because the pressure 
decreases with radius). The inward intrinsic momen- 
tum flux in the absence of rotation (IF„ t < 0) results in 
rotation peaking. The toroidal rotation peaking factor, 
(o/n C)P ) (dftjdr) = (a/fi C)P ) (n mt / Xc ) (l/ mi (R 2 ) s ), as 
a result of this intrinsic momentum transport can be cal- 
culated by assuming II = and i\ = in Eq. ©. 
Fig. 2 shows the toroidal rotation peaking factor with 
Qc,pRa/vti = 0.3 and fl(.ERo/ v u = —0.3. Assuming 
that this dil^/dr is a good estimate for the typical size 
of the rotation gradient, it results in a peak flow in the 
core of up to 40% of the piece of the rotation due to the 
negative radial pressure gradient in the pedestal. 

The existence of two different pinches is one of the 
sources of intrinsic rotation generation according to the 
following restatement of Eq. ([6]) , 
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FIG. 2: The normalized rotation peaking factor generated 
by the intrinsic momentum fluxes in terms of (a) the density 
gradient (R/L n ), (b) the temperature gradient (R/Lt), and 
(c) the safety factor (q). The parameters except the scanned 
variable are the same as in Fig 1. 
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FIG. 3: The normalized profiles of the even parity momentum 
flux : (a) fg™*^ d9-K2 in terms of (v'u/vu), and (b) dv'^n2 
in terms of 9. 



where II mt = U' lnt - ( P < E /^ P ) mi (R 2 } s (n CE - Q c<p ) 
and Il' irlt is the intrinsic toroidal angular momentum flux 
for ri(,_e = 0, f2f, p — 0, and dtt^/dr = 0. 

Result analysis. The toroidal momentum flux can be 
analyzed using the symmetry of the turbulence in 9, v 1 ^ 
and k r (l9| . where k r is the radial wavenumber of the 
turbulence. The momentum flux in the rotating frame 
can be written as IB^ = ^2 k J d9dv'^ir(9 7 v'u , fc r ), and n 
can be expanded in p+ and (B/Bg)p+: ir = tti + ir 2 + ... 
where tt 1 oc f tb ,i(v' ■ C)(V0f • y) and tt 2 oc / t6 . 2 (v' • 
C)(V0i b • y) + / tM (v' • C)(V0 2 b • y) ~ (B/B e )p^ 1 . Be- 
cause the first order fluctuating ion distribution func- 
tion, f t b,i ~ O(/0*/a/), and the corresponding potential 
4>f have the parity f t b,i{0, v'u,k r ) = -f t b,i{-6, -v'u,-k r ) 
and (j) t l'{9,k r ) = —<j) tb (—9 7 —k r ), the toroidal momen- 
tum flux to this order is antisymmetric, Tr{ b (9, v'^ , k r ) = 
— 7r* (— 0, — ufi, — fe r ). As a result, the total sum over tti 
vanishes. 

The different parity of the distribution function 
ftb,2 and the potential 0| & (i.e. ftb,2(8, v'u , k r ) = 
ftb,2(-9i-v[ V -kr) and & b (9,k r ) = ^(-9,-k r )) are 
caused by the Coriolis terms breaking the symmetry in v'^ 
and the additional acceleration term, (ZeRBg/rriic) v' M ■ 
Vr (dftb/de) breaking the symmetry in 9. We fo- 
cus on the additional acceleration term for the case with 
fi^E + il^ yP = because the Coriolis term has already 
been studied [lfj. Using a balance between the accel- 
eration term, (ZeRBg/niic) v' M ■ Vr (df t b,i/de) Q^ iP and 
the time derivative of f t b,2 (dftb,2/dt ~ ftb,il T nh wi th 
r n i the nonlinear decorrelation time [20| ) , we obtain 

ftb,2 ~ — (t; t| +^_ L /2)n c ,p— sm6»-^-, (8) 

where we have used the magnetic drift of concentric cir- 
cular flux surfaces to make the dependence on 9 more 
transparent. The higher order pieces give the even par- 
ity of the next order contribution to the momentum flux: 
7^(0, v'u, k r ) = TT2{—0,—v'n,—k r ). Accordingly, the ratio 



of the momentum flux to the heat flux decreases when the 
safety factor increases because of (Bg/B) ~ (r/R)(l/q) 
(jni is barely modified by the change of the safety factor 
in our simulations). More analysis is needed to com- 
pletely understand the dependence on R/L n and R/ Lt 
in Fig. 2. 

Fig. 3 shows the intrinsic momentum flux as a func- 
tion of v'u and 9 for (Q^pR/vu) — 0.3 and (Q^^eR/vu) — 
—0.3. The dependence on vi and 9 are different from 
those for the Coriolis momentum pinch due to the radial 
electric field by itself. The difference is due to the differ- 
ent mechanisms for the symmetry breaking. The profile 
in v'u in Fig 3-(a) shows that the intrinsic momentum 
transport has a larger contribution from particles with 
the large parallel velocity because f t b,2 has a higher or- 
der moment in v'u due to the magnetic drift effect on the 
acceleration. Physically, the particles with large parallel 
velocity have wider orbits than slow particles and they 
exchange more energy with the radial electric field. The 
asymmetric 9— profile in Fig 3-(b) can be explained by 
the factor sin0 in Eq. (JSJ), which is odd in 9. 

Discussion. The intrinsic momentum flux due to the 
two opposite pieces of the rotation in the pedestal can be 
one of the origins for intrinsic toroidal rotation. When 
ions are not rotating toroidally in the pedestal, a nega- 
tive radial electric field must exist to balance the radial 
pressure drop (£\,e ~ — < 0). The intrinsic mo- 
mentum transport occurs due to the acceleration of the 
particles in the radial electric field that breaks the sym- 
metry of the turbulence. For example, for ATj ~ —IkeV 
in Ar ~ Icm with B e ~ 0.5T, n C:P R a ~ 200km/s. The 
intrinsic momentum transport results in rotation peak- 
ing in the core whose size can be about 40% of f^. p i?0j 
i.e., fl c R ~ 80km/s. 
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